Ferromagnetic planar Josephson junction with 
transparent interfaces: a junction proposal 



D M Heim\ N G Pugach^ ^ M Yu Kupriyanov^, E Goldobin^, 
D Koelle'' and R Kleiner^ 

^ Institut fiir Quantenphysik and Center for Integrated Quantum Science and 
Technology (IQ^"^), Universitat Ulm, D-89069 Ulm, Germany 
^ Skobeltsyn Institute of Nuclear Physics, M. V. Lomonosov Moscow State 
University, 119991 Leninskie Gory, Moscow, Russia 

Faculty of Physics, M. V. Lomonosov Moscow State University, 119991 Leninskie 
Gory, Moscow, Russia 

Physikalisches Institut and Center for Collective Quantum Phenomena in LISA+, 
Universitat Tiibingen, D-72076 Tubingen, Germany 

E-mail: dennis . heim@uni-ulm . de 



Abstract. We calculate the current phase relation of a planar Josephson 
junction with a ferromagnetic weak link located on top of a thin normal metal 
film. Following experimental observations we assume transparent superconductor- 
ferromagnet interfaces. This provides the best interlayer coupling and a low suppression 
of the superconducting correlations penetrating from the superconducting electrodes 
into the ferromagnetic layer. We show that this Josephson junction is a promising 
candidate for an experimental Lp junction realization. 



PACS numbers: 85.25. Cp, 74.78.Fk, 74.45. +c, 74.50.+r 
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1. Introduction 

A ip junction [1, 2] is a Josephson junction with a doubly degenerate ground state, in 
which the Josephson phase takes the values +ip or —ip {0 < ip < tt) [3]. This junction 
being closed into a ring is able to self-generate a fractional flux $o¥'/(2vr), where $o is 
the magnetic flux quantum. 

In this sense the (f junction is a generalisation of the vr junction [4] which has a 
Josephson phase +7r or — vr in its ground state. It has been experimentally demonstrated 
that the it junction improves the performance and simplifies the design of classical and 
quantum circuits [5, 6, 7]. Since the (p junction offers the possibility to choose a special 
value of the phase in the ground state it may further optimize these circuits. 

An experimental evidence of sup junction made of and vr parts [1, 2, 8] was reported 
recently [9]. One half of the junction had the Josephson phase in its ground state and 
the other half the phase it. This was realized [9] by connecting two superconductor- 
insulator-ferromagnet-superconductor (SIFS) junctions in parallel. The advantage of 
this concept is that it is based on the technology already developed for the fabrication 
of O-vr junctions [10, 11]. 

On the other hand this (p junction concept is difficult to realize experimentally 
because, e.g., a step in the thickness of the F layer must be realized with very high 
accuracy [9, 8]. A completely other method, the "ramp-type overlap" (RTO) junction, 
was proposed by Bakurskiy et al. [12]. It only requires one small SFS junction located 
on a thin normal (N) metal layer, see figure 1. This basic setup provides a miniaturized 
ip junction. Moreover, this type of junction has already been realized experimentally for 
the analysis of the double proximity effect [13]. 

A simple model [3] to show that this junction can be used as a junction requires 
its current-phase relation (CPR). By writing it in terms of a sine series 

J(0) = A sin(0) + B sin(20), (1) 
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Figure 1. The geometry of the eonsidered system. The Josephson junction consists 
of two superconducting (S) electrodes separated by a ferromagnetic (F) weak hnk of 
thickness dp and length L. It is located on top of a thin normal (N) metal film of 
thickness d^. 
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where is the Josephson phase, the amphtudes have to obey the conditions [3] 

\B\ > \A\/2 and B < 0. (2) 

The RTO junction, schematically shown in figure 1, can fulfil these conditions 
because the current fiows between the S electrodes through the F metal and the N 
layer. In this way the properties of an SFS and SNS junction are combined. The SFS 
junction can have a negative [14, 15] amplitude Ap in (1), while the SNS junction has 
a positive [14, 16] amplitude A-^ in (1). By adding both the total amplitude A can 
be minimized and a dominant negative amplitude B from the SNS part is obtained to 
fulfil conditions (2). Since supercurrents in SFS junctions are rather small, the SNS 
contribution has to be reduced. This is done by using only a thin normal metal film. 

In the present paper we investigate an RTO junction which has, differently from the 
one proposed in [12], transparent SF interfaces in order to amplify the SFS contribution 
to the total current. This assumption has already successfully been used to describe 
various experiments [17, 18, 13]. As a result, we obtain slightly smaller system sizes 
for the (f junction realization than [12], where weakly transparent interfaces were 
assumed. Moreover, our approach provides a better penetration of the superconducting 
correlations into the F layer which may increase the Josephson current. In the framework 
of transparent SF interfaces we cannot use linearised equations for the SFS part, as it 
was done in [12]. Therefore, we use non-linearised equations in the SFS and SNS part 
for our analytical approach. 

We derive the CPR in the "dirty" limit. For this purpose, we combine the solution 
of the Usadel equations in the N film [12] with the solution of the Usadel equations in 
the SFS layer [19]. The resulting current phase relation consists of three parts: (i) a 
contribution from the SFS layer, (ii) a contribution from the N film and (iii) a composite 
SNFS term. 

The paper is organized as follows. In section 2 we introduce the model of the 
considered Josephson junction in terms of Usadel equations. The analytical expression of 
the CPR of our system is based on this model and presented in section 3. In section 4 we 
use this expression together with realistic system parameters to discuss its applicability 
as ip junction. Finally, an appendix provides a detailed derivation of the composite 
SNFS current. 

2. Model 

The considered Josephson junction is sketched in figure 1. It consists of an SFS junction 
located on a normal metal film. The F layer has a thickness dp and a length L while the 
N layer has a thickness d-^ and is considered as infinitely long. We have chosen the x and 
z axis in directions parallel and perpendicular to the plane of the N film, respectively. 

For the calculation of the current fiowing from one superconducting electrode 
to the other we determine the Green's functions describing our system. Since we 
consider the "dirty" limit, in which the elastic scattering length is much smaller than 
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the characteristic decay length, which allows us to use the Usadel equations [20]. We 
write them in the form [14] 



2 r 

3 



Gj 



dx \ •' dx J dz \ ^ dz 



" = 0, 



G, = ^^^^===, JG{N,F} (3) 

in the N and F layer, respectively. Here, $j and Gj are the Usadel Green's functions 
in the $ parametrization [21]. The frequencies Cj = uj + iH contain the Matsubara 
frequencies u = 7rT(2n + 1) at temperature T, where n = 0, 1, 2, . . ., and the exchange 
field H of the ferromagnetic material which is assumed to be zero in the N layer. The 
decay lengths 



^^ = V2^' ^^=V2^ 
of the superconducting correlations are defined via the critical temperature of the 
superconductor (we use h = = I) and the diffusion coefficients and Dp in the 
normal and ferromagnetic metal, respectively. 

We assume that superconductivity in the S electrodes is not suppressed by the 
neighbouring N and F layers. This assumption is valid in our case of transparent SF 
interfaces with the conditions for the suppression parameters 

-Rbsf^bsf ^ Ps^s ^ , 

7BSF = < 1, 7SF = — ;:- < 1, (5) 

Pf4f Pf4f 

7BSN = T > 7SN = — (6) 

Pn?n Pn?n 

Here, -Rbsn.bsf and Absn.bsf are the resistances and areas of the SN and SF interfaces. 

The values of Pn,f,s describe the resistivity of the N, F, and S metals. 

This allows us to use the rigid boundary conditions [14] 

<l>s(±L/2) = Aexp(±i0/2), Gs = ^=^, (7) 

where A is the absolute value of the order parameter in the superconductor. 
The boundary conditions [21, 22, 14] at the free interfaces are 

^$,=0, JG{N,F}, (8) 
and at the interfaces of the superconductor they are 

7BSNeN^ = ^ [M±m - $n] (9) 



and 



$F = -$s(±V2). (10) 

00 



Additionally we use 

dz Gp V w 



IBNF^F^ = 7^ -$N - $F (11) 
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at the NF interfaces, where 

7BNF = T (12) 

PF^F 

is defined analogous to (6). 

Finally we calculate the total current 

/(0)=/n(0) + /f(0) (13) 

by integrating the standard expressions [14] for the current densities of the N and F 
part over the junction cross section along the z axis. This leads us to 

ttTW 



/n(0) 
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(14) 



x=0 



and 



/f(0) 



ttTW 
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(15) 



x=0 



The width W of the junction along the y axis is supposed to be small compared to the 
Josephson penetration depth. We have chosen the position a; = for the integration 
over the junction cross section since the z component of the current densities vanishes 
there because of the symmetry of the considered junction geometry. 



3. Currents 

In order to calculate the current /(0) from (13) we cannot simply add the current 
through the N layer calculated by Bakurskiy et al. [12] to the SFS current calculated 
by Buzdin et al. [19] because we have to take into account a composite SNFS current 
which appears due to a penetration of superconductivity from the N layer into the F 
layer. Therefore, we split the current If{4>) into a contribution /F,dir(</') due to a direct 
penetration of superconductivity into the F layer and the additional part Jnf(0)- This 
leads us to 

J(0) = /n(0) + /F,dir(0) + /nf(0). (16) 

In the following three sections we derive the expressions of these three currents using 
the scaling 

m = m^- (17) 
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3.1. Current in the N layer 

In this layer we adopt the current 



7 fA\ o^nT r(0) . 



u]>0 

with the definitions 



m = (19) 

27BMf^(v/n^T5V^ + /i(0)) 
5 = — , 7bm = , ^ = -7^, (20) 



7bm 



- c 

-1 



Vuj^ + A^ + l) , (21) 



/i(0) = ^^2+^2^2 cos2(0/2) (22) 

from [12]. Its derivation is based on the assumption L <C and an infinitely long N 
layer. It is calculated with the help of the solution $n(3^) (A9) of the non- linear Usadel 
equations which depends only on the coordinate x because the thickness d-^ <S C,n is 
assumed to be small. 



3.2. Current in the F layer 
The current 

/F,dir(0) = V2 GAdpK e-^^^J^ sin (2kL + ^) sin 0, (23) 

with 

e = — - — ^ ., = v:?Ta^, (25) 

is a result of [19]. It also has been calculated with the help of a solution of the non-linear 
Usadel equations because 7bsf = is assumed. Additionally the condition <^ L is 
required. 



3.3. Composite NF current 

We determine the current /nf(</') by combining the two non-linear solutions $F,dir(a;) and 
$n(3^) of (A6) and (A9) in Appendix A. The main idea is to decompose the ferromagnetic 
Green's function 

$f(x, z) = $F,dir(a;) + $NF(a;, z) (26) 

into a function $F,dir(3;), which corresponds to currents only fiowing in the F layer, and 
a function $NF(a;, z), which corresponds to currents fiowing through the N layer into the 
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F layer. The second function is obtained by linearising the Usadel equations (3) in the 
F layer. Then we connect it to the N layer solution $n(3;) via the boundary conditions. 

The superposition (26) of the solution $F,dir of the non-linear Usadel equation with 
the solution $nf of the linearised Usadel equation is valid because we distinguish in the 
F part between two cases: (i) at x ~ ±L/2 near the boundaries to the S regions the 
Green's function $F,dir is very dominant |$F,dir| ^ I'^'nfI due to a transparent boundary 
between the S and the F part, that is 7bsf = 0; (ii) at a; ~ 0, that is away from the 
boundaries the contribution of $f decays exponentially. Therefore, the contribution 
from the N part is dominant |$nf| ^ |'^'F,dir|- 

As a result (A12) we obtain the current 

16cos(0/2)^F 



Inf{,<P) 



sm 



TBNF^iA^N 

— + -^e"''-^/^ cos (kL 
2 ^/2 V 



x27rT^er(0) sin I 



(27) 



aj>0 



with the definitions of r(0) from (19), k from (24) and O together with rj from (25). 



4. Discussion 



In this section we estimate the geometrical parameters d^^ dp and L, see figure 1, for 
which the considered Josephson junction obeys the ip junction conditions (2). We use 
the analysing scheme of [12] and finally compare our results with the ones obtained 
in [12]. 

We split the sine series amplitudes 

A = A^ + ylF,dir + ^NF, (28) 

B = B^ + Snf (29) 

of the total current (16), scaled according to (17), into parts originating from the current 
of the N layer (18), the F layer (23) and the composite NF current (27). There is no 
amplitude -BF,dir because we have a pure sinusoidal CPR (23) in the F layer. 

In our calculations we chose the temperature T = 0.1 Tc. We make this choice 
because far away from the critical temperature the CPR has larger deviations from the 
sin0 form [16] which results in a larger second harmonic B. As S electrode material we 
chose Nb with = 9.2 K because it is commonly used in superconducting circuits. 

Our first step is to find suitable parameters dp. For this purpose we analyse the 
amplitudes (28) and (29) as a function of L for different values of dp for the same 
parameters as in [12]: d^ = 0.64 ^n, = 0.1 ^at, H = lOTc, A = 1.76 Tc, Pf = Ptss = P 
and 7bnf = 1- Figure 2 shows three typical examples: (a) dp = 0.15 ^n, (b) dp = 0.31 C,n 
and (c) dp = 0.35 ^at. The first (a) and last (c) examples correspond to limiting cases 
where it is difficult to realize a ip junction because the intervals of L where conditions (2) 



Ferromagnetic planar Josephson junction with transparent interfaces 



8 



hold are not large. These intervals of L are highlighted by bold lines. In between the 
two limiting values for dp this line becomes longer. Figure 2 (b) shows an optimum 
situation because there is a wide range of L which yields a (f junction configuration. 




0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 

L L L 



Figure 2. The functions \A\/2 and \B\, based on (28) and (29), as functions of L for 
c^N = 0.64 and three characteristic values of dp- The bold lines correspond to values 
of L where the conditions (2) for the (p junction realization are fulfilled. 



For the optimum value dp = 0.31 we calculate the magnitudes 
Ta = At^WA/{ep) = 0.534 and Tb = Bf^WA/{ep) = -0.106. Inserting them together 
with the amplitude AF,dir from (23) into (2) and neglecting the small NF contributions 
leads us to the condition 



<2\Tb\. (30) 

Here, we use the constant e = ^p/{QAJ^\/hdp) with dp = 0.31 ^n, = 0.0691 and 

^(L) = exp(-2fi:L)sin(2KL + 7r/4). (31) 

From (30) we find the minimum value 0.10 and maximum values 0.17 of L. 

For summarising our suggestion of the geometrical configuration of a <y9 junction we 
use the value .^n = 100 nm for Cu as N layer, a strongly diluted ferromagnet such as FePd 
or the CuNi alloy with = 10 nm and H = 10 Tc as F metal. Our set of parameters then 
become (iN ^ 50 nm, 15 nm ^ dp < 35 nm and 10 nm ^ L < 17 nm, which we compare 
to the values c^n ^ 50 nm, 19 nm ^ dp < 48 nm and 7nm < L < 22 nm of [23]. 

Since we use the same N layer configuration, the value for c^n is the same. But the 
suggested regime for dp differs. A change in this direction was expected because we only 
need a thin F layer since the transparency of our interfaces already amplifies our SFS 
current contribution. The possible range for the length L of the F part is smaller in our 
case but the whole junction configuration is still experimentally feasible. 

5. Conclusion 



We have shown that the considered Josephson junction with a ferromagnetic weak link 
located on a thin normal metal film is a good candidate for a ip junction realization. By 
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choosing transparent SF interfaces we obtained slightly different system sizes for the (p 
junction existence compared to a junction with weakly transparent interfaces. 

The current was split into a contribution through the N layer, the F layer and a 
composite term which described the current flowing through the N and F parts of the 
junction simultaneously. We performed our calculations in the "dirty" limit, that is, the 
currents are obtained from solutions of the non-linear Usadel equations. 

Since our case of a large interface transparency corresponds better to the 
experimental situation [17, 18, 13] than weakly transparent interfaces [23] it is important 
to note that a smaller thickness and length of the F layer have to be chosen than 
predicted in [23]. We are looking forward to experiments realizing this junction and 
its application in classical and quantum devices. 
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Appendix A. NF current derivation 

In this appendix we derive the current (27) which flows through the N and F part of the 
junction, sketched in figure 1, simultaneously. We first linearise the Usadel equations (3) 
and then combine the solution with the Green's functions from [12] and [19]. 

For the linearisation of (3) we assume the superconducting correlation coming from 
the N part into the F part as rather small. Then, the Green's function $nf(x, ^) can 
also be assumed to be small. Using Gnf = sign(a;) we obtain the linearised Usadel 
equation [15] 



* 




(Al) 



with the definitions 



= + i sign(f2)/i, Q 



h 



H 



(A2) 



ttTc 



Its solution in the form of a series 




n=l 



(A3) 



with 




(A4) 



Ferromagnetic planar Josephson junction with transparent interfaces 



10 



and a Fourier coefficient already obeys the boundary condition (8) at the upper 
border {z = + dp). 

The boundary conditions at the left and right end of the F part at x = ±L/2 are also 
already fulfilled. They follow from (10) with 7bsf = 0. Using here the definition (26) of 
$F leads us to the condition 

^pM±L/2) + <I>NF = ^M±L/2). (A5) 
This equation is already fulfilled by the solution 



with 



= 4 arctan 



eexp(±V^^^^ 



(A7) 



e = ^ (A8) 

T] +\u\ + v2?7(?7 + 

from [19] alone. Therefore, the NF Green's function (A3) only has to obey the conditions 
$NF = at a; = ±L/2. Note that we do not need the expression for Gp.dir to finally 
calculate the current. 

In order to obtain the Fourier coefficient 6„ which fixes the solution $nf from (A3) 
we use the boundary condition (11) where we neglect the term <I>nf assuming |$nf| ^ 
|$n|- Now, we replace the Green's function Gp by Gnf = sign(a;) and insert 

$n(x) =rAcos^ + 2i/i(0)r(0)sin( ^ ) ^, = (A9) 



2 ^^"^ ^n' /i(0) 

from [12], where we use the definitions (19), (21) and (22) for r(</)), r and n{4>), 
respectively. By neglecting the real part of $n "we obtain the Fourier coefficient 

2il]L(-l)"r(0)sin(0/2) 
Kn smh(fi;„dF) 4f Tbnft™ 4n 
Our last step is to calculate the current Jnf- Therefore, we insert the Green's 
function (26), which contains the Green's functions from (A3) and (A6), into the 
definition (15) of the F layer current. Due to the condition x = it reduces to a 
sum /f(0) = /nf(0) + -^F,dir(0), whcrc the NF current is defined by 



/nf(0) = i^-- . 



NF<-JF,dir 

02; — 



2epF ,„f-lJdN " 



d d 

$F,dir(a;) -^^pi-Uj) - $^,dir(-^) ^'^Nf(w) 



(All) 



x=0 

and /F,dir(0) is the current flowing only through the F layer [19] summarized in (23). 

We insert the Green's functions $nf and $F,dir from (A3) and (A6) into (All). 
Then, by finally using the the approximation Q ^ ih, which holds for the condition 
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ttTc <^ H, we obtain the scaled current 




cos { kL + — 
V 4 



kL/2 



) 



(A12) 
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